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Abstract 

Background: The mechanics of the right ventricle (RV) are not well understood as studies of the RV have been 
limited. This is, in part, due to the RV's thin wall, asymmetric geometry and irregular motion. However, the RV plays 
an important role in cardiovascular function. This study aims to describe the complex mechanics of the healthy RV 
using three dimensional (3D) cine displacement encoding with stimulated echoes (DENSE) cardiovascular magnetic 
resonance (CMR). 

Methods: Whole heart 3D cine DENSE data were acquired from five healthy volunteers. Tailored post-processing 
algorithms for RV mid-wall tissue tracking and strain estimation are presented. A method for sub-dividing the RV 
into four regions according to anatomical land marks is proposed, and the temporal evolution of strain was 
assessed in these regions. 

Results: The 3D cine DENSE tissue tracking methods successfully capture the motion and deformation of the RV at 
a high spatial resolution in all volunteers. The regional Lagrangian peak surface strain and time to peak values 
correspond with previous studies using myocardial tagging, DENSE and strain encoded CMR. The inflow region 
consistently displays lower peak strains than the apical and outflow regions, and the time to peak strains suggest 
RV mechanical activation in the following order: inflow, outflow, mid, then apex. 

Conclusions: Model-free techniques have been developed to study the myocardial mechanics of the RV at a high 
spatial resolution using 3D cine DENSE CMR. The consistency of the regional RV strain patterns across healthy 
subjects is encouraging and the techniques may have clinical utility in assessing disrupted RV mechanics in the 
diseased heart. 
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Background 

Right ventricular (RV) function may be impaired in a 
number of heart conditions, including myocardial infarc- 
tion, congenital heart disease and cardiomyopathy [1]. 
The function of the RV may also be affected in diseases 
of the left ventricle (LV) where it is difficult to ignore 
the complex nature of ventricular interaction [2,3]. In 
the past, the importance of the LV in cardiac research 
has overshadowed the study of the RV. This neglect is, 
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in part, because the RV is difficult to image. The wall of 
the RV myocardium is thin (2-5 mm) when compared 
to that of the LV (7-11 mm) [2]. Furthermore, the RV 
has a complex geometry, eccentric motion [1,4] and it is 
heavily trabeculated, thus it does not offer the clearly 
defined endocardial margins typically seen in the LV. 

Various cardiovascular magnetic resonance (CMR) 
studies have assessed properties such as strain, motion 
and volumes of the RV. The CMR techniques used 
include turbo gradient echo and steady state free preces- 
sion (SSFP) [5], myocardial tagging [4,6,7], phase con- 
trast velocity encoding [8], strain encoded (SENC) CMR 
[9] and 2D multi-slice displacement encoding with sti- 
mulated echoes (DENSE) [10]. 
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The standard cardiac imaging planes are based on the 
relatively uncomplicated geometry of the LV. However, 
no standardized localization exists for the RV. Two 
dimensional strain estimates of the LV are reliable in 
short axis views because the through plane motion is 
relatively uniform [11]. This does not apply to the RV 
because of its complex shape and motion, so 2D ima- 
ging may be insufficient to accurately assess strain in 
the RV [4], To date, no detailed studies of the RV have 
been presented using DENSE CMR. 

DENSE is a quantitative CMR technique used for 
measuring myocardial displacement and strain. DENSE 
encodes tissue displacement directly into the image 
phase (typically with reference to end diastole), thus 
allowing for the extraction of motion data at a pixel 
resolution [12,13]. A recently developed free breathing 
3D cine DENSE sequence [14] is well suited for quanti- 
fying the complex behavior of the RV. A high signal to 
noise ratio (SNR), which is necessary for imaging the 
thin RV walls, is achieved using a spiral k-space trajec- 
tory and three point phase cycling. Furthermore, 
DENSE is inherently a black-blood technique, which 
together with fat suppression, provides a more accurate 
delineation of the RV walls. 

The purpose of this study was to develop tailored pro- 
cessing techniques for assessing detailed 3D RV motion 
and surface strain using 3D cine DENSE, and to quan- 
tify these parameters for the healthy human heart. 

Methods 

Imaging protocol 

Whole heart 3D cine DENSE data were acquired from 
five healthy male volunteers (age range 21 - 45) on a 
1.5T MRI scanner (Siemens MAGNETOM Avanto, 
Erlangen, Germany) using a two-channel anterior body 
array coil and an eight-channel spine array coil. 

The entire heart was imaged with the imaging volume 
aligned along the heart's long axis and at a 2.8 x 2.8 x 5 
mm 3 spatial resolution and 32 ms temporal resolution. 
Displacement was encoded in three orthogonal direc- 
tions and a spiral k-space trajectory was used with 6 
interleaves per 3D partition. Spiral aliasing artifacts 
existed on the edge of images, but not on the heart, 
which was placed at the center of the field of view. 
Images were acquired during a 20.5 ± 5.7 min scan time 
with prospective ECG gating, diaphragmatic navigator 
respiratory gating, and no multichannel parallel imaging. 
Other imaging parameters include: field of view (FOV) 
= 360 x 360 x 140 mm 3 (in-plane dimensions defined 
by the diameter of the spiral), displacement encoding 
frequency k e = 0.06 cycles/mm, ramped flip angle up to 
20 degrees, TR = 16 ms, TE = 1.3 ms, cardiac phases = 
22 and number of acquired 3D partitions = 14 (zero- 
padded to 28 during image reconstruction). After the 



Fourier Transform in the partition direction, three parti- 
tions at each end of the volume were discarded to avoid 
aliasing. The resulting 3D image matrix size was 128 x 
128 x 22. The number of partitions chosen per subject 
in order to include the RV was approximately 14, which 
corresponds to an overall stack thickness of 70 mm. All 
imaging was conducted after informed consent and in 
accordance with protocols approved by the University of 
Virginia Institutional Review Board. 

Post processing 

A number of tailored post processing steps were imple- 
mented to analyze the RV myocardial DENSE data. All 
software development was performed using MATLAB 
(The Mathworks, Natick, MA). 

Contouring, phase unwrapping and tissue tracking 

Both the RV and LV myocardium was demarcated from 
surrounding structures by manually drawing each epi- 
cardial and endocardial contour on the reconstructed 
short-axis DENSE magnitude images. 

Spatio-temporal phase unwrapping [15] was per- 
formed within the contoured area to remove phase 
aliasing, and the phase data from the three encoding 
directions were combined to create 3D Eulerian displa- 
cement fields. Figure 1 illustrates contoured cine 
DENSE magnitude and phase-unwrapped images and 
the corresponding 3D displacement fields at end systole 
for a single short-axis slice. 

There are typically too few transmural pixels span- 
ning the RV to acquire a full 3D strain tensor, so 2D 
strain at the RV mid-wall surface was calculated. The 
mid-line between epicardial and endocardial contours 
was used to produce tissue tracking seed points, which 
were spaced at pixel-distance intervals on the first car- 
diac phase. The mid-line contours were spatially 
smoothed within and across partitions using 4 th order 
polynomial fitting to ensure a continuous RV surface. 
An example of the RV mid-line points is shown in 
green in Figure 1(a). 

A 2D tissue tracking algorithm involving interpolation 
of the displacement fields [15] was directly extended to 
3D for this application. The position of each mid-line 
point along its motion trajectory throughout the cardiac 
cycle was estimated using 3D distance weighted linear 
interpolation. The full 3D mid-line motion trajectories 
were calculated by subtracting the interpolated vectors 
of successive frames from one another. Prospective gat- 
ing only allowed about 90 percent of the cardiac cycle 
to be imaged, therefore, temporal fitting was done for 
each ordinate direction of each of the trajectories using 
a 10 th order polynomial, as 5 th order Fourier basis func- 
tions have been shown to sufficiently describe cardiac 
motion [16]. 



Auger et al. Journal of Cardiovascular Magnetic Resonance 2012, 14:4 
http://www.jcmr-online.eom/content/14/1/4 



Page 3 of 9 




Figure 1 A single short axis slice of the 3D cine DENSE data at mid systole, (a) DENSE magnitude image illustrating the myocardial mid- 
line for the RV (green dotted line), (b, c, d) DENSE unwrapped phase images encoded for motion in the x, y and z directions, respectively, (e) 
Corresponding 3D DENSE displacement field. Red lines correspond to epicardial surfaces and blue lines correspond to endocardial surfaces. 



Lagrangian surface strain 

Lagrangian strain was calculated directly from the 3D 
motion trajectories and oriented tangential to the RV 
mid-line surface. Considering an arbitrary point on the 
mid-wall, the deformation gradient tensor F was calcu- 
lated using two to four of the nearest neighboring mid- 
line points as follows 



F = [V / 1 V^]/[ViV 2 V 3 



(1) 



where Vi, V2 and V 3 represent the pre-deformation 
vectors which are perpendicular to each other in 3D 
space as shown in Figure 2 and form the columns of 
one matrix. V\, V 2 and V' 3 represent the corresponding 
post-deformation vectors and form the columns of 
another matrix, and the "/" operator is a matrix right 
division. 

In practice, depending on the position of the mid-line 
point on the RV wall, each motion trajectory could have 
adjacent points to the right and/or left, and on 



partitions above and/or below. The deformation gradient 
tensor is therefore calculated for a minimum of two and 
a maximum of four co-planar vector configurations for 
each mid-line point. An example of four co-planar vec- 
tor configurations is illustrated in Figure 2a, and an 
example of two co-planar vector configurations pre- and 
post-deformation is illustrated in Figure 2b and 2c, 
respectively. Orthogonal unit vectors (e.g. V3 and V'3) 
were created perpendicular to the RV surface using 
cross products of each set of co-planar vectors (e.g. VI 
and V2, V'l and V'2). Three dimensional Lagrangian 
strain was then calculated for each deformation gradient 
configuration as follows 



1 



(F T F — I) 



(2) 



where T denotes the transpose operation, and I repre- 
sent the identity matrix. Note that this construction of a 
3D tensor using a unit normal vector assumes that the 




Figure 2 2D Lagrangian surface strain estimation, (a) RV surface image with superimposed mid-line points (white dots) representing the 
partitions. All possible vector configurations are shown for a chosen point (C) using mid-line points (A)bove, (B)elow, to the (L)eft and to the (R) 
ight. (b, c) Pre- and post-deformation, respectively, for a single vector configuration. Strain for the chosen point C is calculated using the 
configuration of points above and to the right (ACR). and V 2 represent the co-planar vectors formed. V 3 represents the orthogonal unit 
vector calculated from the cross product of and V 2 . and V 2 represent the corresponding post deformation co-planar vectors, and V' 3 is 
the cross product of V' n and V 2 . 
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cardiac tissue is non-compressible and has been imple- 
mented here to simplify the strain calculations. The 3D 
strain tensor was decomposed into its corresponding 
eigenvalues (strain) and eigenvectors (direction of defor- 
mation). The eigenvalue corresponding to the orthogo- 
nal unit vector (E]J is aligned with the direction normal 
to the surface (radial) and is zero because of the afore- 
mentioned unit normal vector assumption. The remain- 
ing two eigenvalues (E 2 and E 3 ) provide an estimate of 
the 2D tangential surface strain. The eigenvalues E 2 and 
E 3 are of similar magnitude, and there are regional dif- 
ferences in the myofiber (and eigenvector) orientation 
across the RV surface, so we were unable to discern 
between E 2 and E 3 . The results of the surface strain ana- 
lysis are thus presented as a mean Lagrangian 2D princi- 
pal strain, which has been averaged for the various co- 
planar vector configurations. 

It has been shown that tangential strain in the circum- 
ferential and longitudinal directions are good indicators 
of LV dyssynchrony [17]. In order to infer details about 
the strain directions in the RV, one dimensional (ID) 
strains were calculated in the longitudinal and circum- 
ferential directions in the same manner as usually done 
for the LV, and relative to the orientation of the LV. 
Each ID strain was calculated from the motion trajec- 
tories by only considering points above and below (long- 
itudinal) or to the left and right (circumferential) of each 
mid-line point. This is also illustrated in Figure 2. 

Anatomical sub-divisions for strain-time analyses 

Unlike the LV, there is currently no standard method 
for dividing the RV into anatomical sub-regions. Pre- 
vious studies of the RV using myocardial tagging have 
adopted different conventions for sub-dividing the RV. 
Klein et al. divides each RV basal, mid and apical short 
axis slice into three regions: superior, mid and inferior 
[18]. Haber et al. [4], divided the RV into four regions 
according to anatomical landmarks. The parietal and 
septal bands were used to demarcate the outflow region, 
while the free wall was identified based on the normal- 
ized height of the septum in a long axis view [4]. Fayad 
et al. divides the RV using the supraventricular crest 
(SC), moderator band, papillary muscles and the tricus- 
pid valve as landmarks, creating four distinct regions: 
the inflow, outflow, mid and apical regions [19]. 

We propose dividing the RV surface in a similar man- 
ner to Fayad et al. These regions are illustrated in Figure 
3. The divisions were defined by navigating through the 
3D DENSE magnitude-reconstructed images and identi- 
fying three landmark coordinates as follows: 

Coordinate 1: The SC is an accentuation of the mus- 
cular wall demarcating the inflow and outflow tracts 
within the RV. Navigate from apex to base along the 
short axis magnitude images until the SC appears in the 



image. The first coordinate is marked on the RV free 
wall at the level of and adjacent to the SC. 

Coordinates 2 and 3: The moderator band is a muscu- 
lar band connecting the interventricular septum to the 
anterior papillary muscle of the tricuspid valve. Navigate 
from base to apex along the short axis magnitude 
images and define the slice where the moderator band 
becomes evident as a segment of myocardium crossing 
the RV cavity. The second and third coordinates are 
defined as the anterior and inferior RV-LV insertion 
points at this level, respectively. 

The apical RV region is defined as the myocardium 
lying apical to the slice identified by the moderator 
band. The inflow, mid and outflow RV regions are 
defined using two planes: (i) a plane defined by coordi- 
nate 1, coordinate 2 and a point displaced normal to the 
RV mid-line surface in this region, and (ii) a plane 
defined by coordinate 1, coordinate 3 and a point dis- 
placed normal to the RV mid-line surface in this region. 
The strain values (2D and ID Lagrangian) were aver- 
aged in each of the 4 segments and assessed through 
the duration of the cardiac cycle for all 5 volunteers. 

Results 

The frame-to-frame 3D mid-line motion trajectories for 
both the LV and RV are illustrated in Figure 4. The red 
and blue colors represent the positions along the motion 
trajectories for the RV and LV, respectively. The evolu- 
tion of displacement during the cardiac cycle can be 
appreciated in detail, and the large range of motion of 
the RV compared to the LV is evident. 

The regional mean RV midwall strain-time curves for 
the 5 volunteers are plotted in Figures 5 and 6. Figure 5 
shows the mean of the principal strains E 2 and E 3 , while 
the ID tangential circumferential and longitudinal 
strains are shown in Figure 6. The data are represented 
as mean ± one standard deviation. A summary of the 
maximum strain and the time to peak strain for each 
RV anatomical region is given in Table 1. 

All strain values are consistently negative, indicating 
muscle shortening or contraction, and the longitudinal 
shortening is greater for all regions when compared to 
circumferential shortening. 

The inflow region has the lowest peak strain value of 
all the regions, while the outflow region generally 
reaches the highest peak strain. However, Figure 6 illus- 
trates that during tangential shortening in the circum- 
ferential direction, the apical region tends to have a 
slightly greater strain than the outflow region. Good 
consistency between the patterns of strain evolution is 
shown for the various strain estimates in Figures 5 and 
6. The four regions each arrive at a maximum strain at 
slightly different time points in the cardiac cycle. The 
time sequence of regional myocardial shortening, as 
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(c) 



Figure 4 3D mid-line motion trajectories, (a, b) show the RV (red) and LV (blue) positions during end-diastole and end-systole, respectively, (c, 
d and e) represent the motion trajectories for the respective slices 1,2 and 3. 
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Figure 5 Mean principal 2D Lagrangian regional RV surface strain-time curves for 5 normal volunteers. 



confirmed by all strain estimates is: inflow, outflow, mid 
then apical 

Figure 7 demonstrates 3D cine DENSE bi-ventricular 
midwall tissue tracking in two views for one volunteer 
at end-diastole, mid-systole and end-systole, with the 
tracked points colour coded according the mean princi- 
pal strain. The additional files 1 and 2 show these data 
for all cardiac phases in plane and oblique views, 
respectively. 

Discussion 

This paper introduces post processing algorithms devel- 
oped for 3D tissue tracking and strain analysis of the RV 
from 3D cine DENSE data. Unlike previous studies of 3D 
RV motion using tagging [4,7], no models were used to 
describe the underlying data. Furthermore, no spatial 
smoothing was applied to the data, and the only form of 
filtering was the temporal polynomial fitting applied to 
the mid-line motion trajectories. This temporal fitting 
provides a more realistic trajectory behavior and allows 
for motion to be estimated at any temporal resolution. 

Both the peak strain and the time to peak strain were 
found to vary across the different RV regions. The 
inflow region consistently demonstrated lower strains 
and early strain peak times, while the apex consistently 
showed larger strain values and later peak times. Some 
variation does exist in the outflow region when compar- 
ing the different strain calculation types. 



As there is no current standard for dividing the RV, 
and multiple approaches for calculating strain have been 
adopted, it is difficult to make a direct comparison with 
results from previous studies. In general, however, the 
magnitude trends and timing parameters of the strain 
measurements from this study compare favorably with 
previous RV studies using myocardial tagging, SENC 
and DENSE. From our study the peak mean principal 
strains values in the apical, mid and outflow regions are 
-0.17, -0.16 and -0.18, respectively. These values are 
similar to the corresponding values of -0.17, -0.19 and 
-0.22 reported in a tagging study by Haber et al. [4]. 

In the SENC study by Hamdan et al. [9], they show 
longitudinal strains to be largest at the base, with a 
lower strain at the apex, and the lowest strain in the 
mid region. For circumferential strain they show the lar- 
gest strain at the apex, with a lower strain in the mid 
and basal regions. These regional strain comparisons 
show excellent agreement with those presented in Table 
1. In our results, during longitudinal contraction, the 
largest strain is found in the outflow region (basal), fol- 
lowed by the apical and mid regions, while circumferen- 
tial strain in the apical region is larger, followed by the 
mid and outflow (basal) regions. 

The ID peak longitudinal strains in Table 1 are all lar- 
ger than the corresponding ID peak circumferential and 
mean principal strains, a result confirmed in tagging 
studies by both Haber et al [4] and Fayad et al. [19]. 
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Figure 6 1 D Lagrangian regional RV surface strain-time curves for 5 normal volunteers. 



Table 1 Right ventricular regional Lagrangian strain and 
time to peak strain 



Strain quantification method and 
anatomical region 


Maximum 
strain 


Time to 
peak (ms) 


1 D Circumferential 






Inflow 


-0.10 ± 0.04 


224 


Outflow 


-0.15 ± 0.05 


256 


Apex 


-0.17 ± 0.06 


320 


Mid Region 


-0.16 ± 0.06 


320 


ID Longitudinal 






Inflow 


-0.16 ± 0.06 


256 


Outflow 


-0.22 ± 0.08 


288 


Apex 


-0.20 ± 0.07 


352 


Mid Region 


-0.18 ± 0.06 


288 


Mean principal strain 






Inflow 


-0.12 ± 0.05 


224 


Outflow 


-0.18 ± 0.07 


256 


Apex 


-0.17 ± 0.06 


320 


Mid Region 


-0.16 ± 0.05 


288 



Fayad et aL, who used a very similar method to ours for 
sub-dividing the RV, found regional peak circumferential 
strains to range from largest to smallest in the following 
order: apical, mid, inflow and outflow. For peak longitu- 
dinal strains the corresponding order was: outflow, 
inflow, apical and mid regions. The peak strain values in 
Table 1 also correlate well with the segmental shorten- 
ing presented by Fayad et al., except for the inflow 
region where our values are consistently lower. 

Wen et al. [10] present circumferential strain results 
from 2D DENSE data as an average of the entire RV 
free wall, in three slices: basal, mid and apical. Here, the 
strain is derived from the lengthening of the RV con- 
tour. They also found decreasing peak circumferential 
strain from apex to base, but the results plotted in [10] 
have a consistently larger magnitude than those pre- 
sented here. 

Limitations of our study include a lengthy scan time, 
the need for manual contouring, and the difficulty to 
discern between the minimum and intermediate 



Auger et al. Journal of Cardiovascular Magnetic Resonance 2012, 14:4 
http://www.jcmr-online.eom/content/14/1/4 



Page 8 of 9 




Figure 7 Tissue tracking and mid-wall strain. LV and RV mid-line images representing mean principal 2D Lagrangian surface strain in two 
views, (a, d) End-diastole, (b, e) mid-systole and (c, f) end-systole. The dots represent the motion trajectory positions, while color represents 
mean principal strain. 

\ J 



principal strains. A reduction in the scan time could be 
accomplished by incorporating parallel imaging and/or 
outer volume suppression. The time taken to contour 
the images may be considerably reduced with the exten- 
sion of previous work on motion guided segmentation 
for 2D cine DENSE images [20]. Even with the use of 
the eigenvectors, discerning between the two principal 
strains was not possible. A geometric or finite element 
model of the RV could solve this problem, whereby the 
equivalent of an RLC coordinate system for the LV is 
defined according to the RV geometry. 

The 2D surface strain estimates are expected to have a 
lower accuracy in the through-plane direction because 
the through-plane image resolution is lower than the in- 
plane resolution. This should lead to a less accurate 
representation of strain in the longitudinal direction 
compared to the circumferential direction, but this is 
not evident in the standard deviation curves in Figure 6. 

An inherent limitation of DENSE is the reduction in 
SNR associated with the stimulated echo. The SNR was 
typically sufficient to clearly identify the RV because of 
the short echo times associated with the spiral readout. 
However, signal loss due to through-plane dephasing is 
particularly pronounced in the RV because of its large 
through-plane strain [15]. This reduction in SNR was 
evident in a few instances, and these data had to be 
excluded from the analyses. For future work on the RV, 
this through-plane dephasing could be minimized by 
lowering the displacement encoding frequency, but at 
the expense of displacement sensitivity. 



Conclusions 

Three dimensional cine DENSE CMR is a promising tech- 
nique for assessing RV mechanics. The methods presented 
here enable the motion and deformation of the entire RV 
to be captured at a pixel resolution. The RV mechanics 
measured by 3D cine DENSE is comparable to that of pre- 
vious studies using tagging, SENC and 2D DENSE. 

We further describe a method for dividing RV into 
regions according to anatomical land marks. The differ- 
ing temporal strain evolution in the proposed RV anato- 
mical sub-divisions indicates that these regions may also 
be functionally distinct. Future work will involve apply- 
ing these techniques to studying variations in regional 
RV function in cardiac disease. 

Additional material 



Additional file 1: Three dimensional mid-line motion trajectories in 
a plane view of the LV and RV seeded from the mid-line points, for 
each partition/slice. Motion trajectories begin at end diastole. The dots 
represent the motion trajectory positions, while the color represents 
mean principal 2D Lagrangian surface strain. 

Additional file 2: Three dimensional mid-line motion trajectories in 
an oblique view of the LV and RV seeded from the mid-line points, 
for each partition/slice. Motion trajectories begin at end diastole. The 
dots represent the motion trajectory positions, while the color represents 
mean principal 2D Lagrangian surface strain. 
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